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Using the holographic gauge-gravity duality, we find a solution for an isolated vortex and a vortex 
lattice in a 2+1 dimensional p-wave superconductor, which is described by the boundary theory 
dual to an SU(2) gauge theory in 3+1 dimensional anti-de Sitter space. Both p x + ip y and p x — ip y 
components of the superconducting order parameter, as well as the effects of a magnetic field on these 
components, are considered. The isolated vortex solution is studied, and it is found that the two 
order parameter components have different amplitudes due to the time reversal symmetry breaking. 
The vortex lattice for large magnetic fields is also studied, where it is argued that only one order 
parameter component will be nonzero sufficiently close to the upper critical field. The upper critical 
field exhibits a characteristic upward curvature, reflecting the effects of field-induced correlations 
captured by the holographic theory. The free energy is calculated perturbatively in this region of the 
phase diagram, and it is shown that the triangular vortex lattice is the thermodynamically preferred 
solution. 



Recently the gauge-gravity duality [TJ[2] nas provided a 
new means by which to explore strongly interacting the- 
ories in condensed matter systems. One of the earliest 
and most successful of such models has been shown to 
exhibit the key properties of superconductivity: a phase 
transition at a critical temperature, where a spontaneous 
symmetry breaking of a U(l) gauge symmetry in the 
bulk gravitational theory corresponds to a broken global 
1/(1) symmetry on the boundary, and the formation of 
a charged condensate. Gravity duals have so far been 
found for s-wave superconductors [3J 2| (in which the 
Cooper electron pairs have angular momentum I = 0) , as 
well as for p-wave (I = 1) [Ml] and d-wave (1 = 2) [HI E] 
superconductors. The field of holographic superconduc- 
tivity has since grown rapidly (see Ref. |10j for reviews). 
Such a dual description provides a window through which 
we might hope to obtain insight into the properties and 
behaviors of superconductors and superfluids that defy 
description by more traditional approaches. 

The basic recipe for creating a holographic supercon- 
ductor involves the introduction of a black hole and a 
charged scalar field into anti-de Sitter (AdS) spacetime 
in d+1 dimensions. According to the AdS/CFT corre- 
spondence, this theory is dual to a d-dimensional field 
theory that exists on the boundary of this space, where 
the boundary value of the bulk field is related to the 
expectation value of an operator, which is interpreted 
as the superconducting order parameter, in the bound- 
ary theory, and the temperature of the boundary theory 
is given by the Hawking temperature of the black hole. 
Below a critical temperature T c the field condenses, and 
the operator on the boundary acquires a nonzero expecta- 
tion value, which corresponds to a superconducting phase 
transition. 

The simplest case in the above scenario is that in which 
the electromagnetic and scalar fields are a function only 
of the radial AdS coordinate, having a boundary value 
which is spatially uniform. Such solutions were soon ex- 
tended to include solutions for isolated vortices, which 



feature a spatially non-uniform order parameter, in s- 
wave superconductors [TT1 - [T3] and vortex lattices in s- 
wave [21] and <i-wave [14] superconductors. Extending 
these results to p-wave superconductors is of interest for 
a number of reasons. The order parameter in these sys- 
tems has multiple components and breaks time-reversal 
symmetry, which leads to a richer set of possibilities than 
is possible in the simpler s-wave superconductors. In ad- 
dition, the fact that the gravity dual is a SU(2) gauge 
theory makes its development more straightforward than 
for G?-wave superconductors, where the dual theory in- 
volves a spin-2 field in a gravitational background. The 
theory of holographic p-wave superconductors is also at- 
tractive because it has fewer free parameters than the 
charged scalar field theories that describe s- and d-wave 
superconductors. 

The action proposed in Ref. [7] to describe a p-wave 
holographic superconductor is 
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where k is the gravitational coupling, R is the Ricci scalar 
curvature, L is the radius of AdS space and F* = d^A^ — 
d v A° + e^A^Al is the SU{2) Yang-Mills field strength. 
The bulk gravitational theory is described by the AdS- 
Schwarzschild metric: 
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where z is the radial AdS coordinate and g(z) = 1 — 
(z/zo) 3 . The black hole horizon at z = z is related 
to the Hawking temperature of the black hole, which is 
equal to the temperature of the boundary theory, by zq = 
3/(47rT). Our calculations will be performed in the probe 
limit, in which there is no back-reaction of the gauge field 
on the metric [3]. 

As a starting point we consider the following ansatz 
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for the gauge field: 



energy solution, yields 
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Here r a are the generators of SU(2), which obey the 
relation [r a , r b ] — e abc T c and are related to the Pauli 
matrices by r a = a a /2i. Following Refs. [5]-[7], we inter- 
pret the U(l) subgroup generated by r 3 as the group of 
electromagnetism, so that &(x, y, z) and A 3 y (x, y, z) are 
the electromagnetic scalar and vector potentials, respec- 
tively. Because r 3 generates a rotation in the 1-2 "plane," 
which is the analogue of a rotation in the complex plane 
in an ordinary Ginzburg-Landau theory of superconduc- 
tivity, the (real) scalar fields w± (x, y, z) are charged un- 
der this U(l), and they represent the amplitudes of the 
p x ± ip y components of the superconducting order pa- 
rameter, respectively. 

To study the case of an isolated vortex, we switch 
boundary coordinates from {x,y) to (r, ((>). In an or- 
dinary Ginzburg-Landau theory with a complex order 
parameter, the vortex solution is found by replacing 
ip(r,<j)) — ► e m ^tp(r), such that the phase changes by 2im 
as one goes around the vortex core, and n is known as 
the winding number. By analogy, the vortex ansatz for 
the p x ± ip y superconductor is given by the replacements 
w±(r, (f>, z) — ¥ exp(2n±(j)T 3 )w±(r, z) in Eq. ([3]), where n± 
are the (integer) winding numbers for the two compo- 
nents of the superconducting order parameter. With this 
modification, the gauge field ansatz becomes 



A x = w+(r, z) cos(n + 0) + w_ (r, z) cos(n_( 
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A 2 , = w+(r, z) sm(n+<f>) — W-(r, z) sin(n_ 
A 2 = w + (r, z) cos(n + <p) — w_ (r, z) cos(n_ 
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Furthermore, we assume that the electromagnetic scalar 
and vector potentials are rotationally symmetric and 
given by 3>(r, z) and A?(r, z), respectively. 

The next step is to determine and numerically solve 
the equations of motion for this ansatz. The Yang-Mills 
equations are 
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To get the equation of motion for w + (r, z) , we add the 
Yang-Mills equations ^ for {a,v) = (1, as) and (a,i>) = 
(2,y). We find that to obtain consistent equations, with 
w + independent of (j>, requires that the winding numbers 
are related by n_ = n + + 2. This relation between the 
winding numbers is also seen in the ordinary Ginzburg- 
Landau theory of p-wave superconductors [16] , and is ul- 
timately due to the presence of mixed gradient terms such 
as (D x + iD y ) 2 in the full equations of motion. Choos- 
ing n± = =pl, which we expect corresponds to the lowest 
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Similarly, the equations of motion for w- and the gauge 
fields are 
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=d z {gd z A%) + rd r [~d r A?A + rd r (w 2 _ - w\) 
+ (w + + W-)d r Al + r(w + + w_)(w 2 _ - w 2 + ). 
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These equations can be solved numerically, subject to 
appropriate boundary conditions. At the boundary of 
AdS at z = 0, the fields have the limiting forms 



w± ={0±)z- 
$ =/i — pz -f 

B(r) =-d r A% 
r v 



(10) 



where {0±) are interpreted as the two components of 
the superconducting order parameter (the square root is 
necessary since 0± has mass dimension 2, whereas the 
superconducting order parameter should have mass di- 
mension 1), p is the chemical potential, p is the charge 
density, and B{r) is the magnetic field. In our numerical 
solution, we specify the value of p, as well as the con- 
ditions w±(z = 0) = and d z A^ — 0. As discussed in 
Ref. |13| , this last condition is a Neumann boundary con- 
dition, which, unlike the more commonly used Dirichlet 
condition, allows for the presence of a dynamical gauge 
field in the boundary theory, and also, according to the 
AdS/CFT dictionary, leads to vanishing of the current 
operator in the boundary theory. At the horizon (z = Zq), 
we require $ = and A 3 ^ be regular. At the vortex core 
(r = 0), the boundary conditions are w± — 0, A 3 ^ — 
and 9, $ = 0. Finally, far from the vortex core at the edge 
of our solution domain (r — R), we require d r w± = 0, 
d r <& = and A 3 ^ — 1. This last condition ensures that 
there is one quantum of magnetic flux passing through 
the vortex region (15) . with Jd 2 rB(r) = 2ir. To obtain 
our numerical solutions we have used the COMSOL 3.4 
package [T7] . 
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FIG. 1. (Top) Spatial profile of the w+ (solid) and W- 
(dashed) components of the superconducting order param- 
eter for an isolated vortex in a magnetic field at temperature 
T/fj, — 0.032. (Bottom) Magnetic field profile for the same 
vortex configuration. 



Fig. [T] shows the spatial profile of the two components 
of the superconducting order parameter for an isolated 
vortex, through which a single quantum of magnetic flux 
penetrates the superconductor. It can be seen that, as 
the order parameter approaches its bulk value far from 
the vortex core, the w_ component has a slightly smaller 
amplitude than the w+ component. This is a conse- 
quence of the breaking of time-reversal symmetry, and 
the fact that the two components do not couple to the 
external field in the same way. The difference between 
the two components will grow as the field increases, how- 
ever vortices will tend to proliferate at higher fields, so 
that the picture of an isolated vortex eventually ceases to 
be valid. The hump at r ~ 0.2 is an interesting feature 
that appears to be present at all temperatures. It may 
arise from the fact that the gradient terms in our theory 
are different from those in the usual Ginzburg-Landau 
theory, where such a hump is generally not present [16] . 

The lower part of Fig. [T] shows the profile of the mag- 
netic field near the vortex core. The exponential decay 
of the field with distance from the vortex core is a gen- 
eral property of superconductors and is also seen in the 
conventional Ginzburg-Landau theory. It is interesting 
to compare the size of the vortex core (the so-called "co- 
herence length") £ ~ 0.1 to the penetration depth of the 
magnetic field A ~ 3 [18]. The ratio of these quantities 
defines the Ginzburg-Landau parameter k = A/£ ~ 30. 



The fact that k ^> 1 means that the holographic p-wave 
superconductors are strongly type II, and are therefore 
expected to exhibit a vortex lattice solution near an up- 
per critical magnetic field H C 2 . Most of the superconduc- 
tors that attract widespread theoretical interest, includ- 
ing the high-temperature cuprate superconductors, fall 
in this regime. 

We begin our discussion of the superconducting prop- 
erties near the upper critical field by considering the case 
of an ordinary, non-holographic p-wave superconductor 
in a magnetic field. In this case the ground state that 
is realized is known to depend on the coefficients of the 
kinetic terms in the free energy. For a two-dimensional 
superconductor with two complex components, the ki- 
netic terms allowed by symmetry are Ki(Dir)j)* Dirjj, 
K-2(Dirii)* DjT)j and K 3 (Diijj)* Djrji, and the ground 
state that is realized depends on the values of these coef- 
ficients [IS]. For (K2 + K 3 )/Ki < the free energy is not 
bounded from below and the theory is unstable. (This 
is true at least for the linearized version of the theory. 
The theory may still be stable when higher order terms 
are included.) In the stable region, one of two possible 
ground states is realized. For (K2 + K$)/K\ > and 
K 2 -K 3 > {K 2 + K 3 ) 2 /(2K 1 +K 2 + K 3 ) only one or- 
der parameter component is nonzero (e.g. rj + ~ rj\ + ir/2 
is nonzero if the field is in the +z direction), and this 
component is in the lowest (n = 0) Landau level. Such 
a state shall be denoted as |0) + . On the other hand, for 
K 2 -K 3 < (K 2 + K 3 ) 2 /{2K 1 +K 2 + K 3 ), both order 
parameter components are nonzero, and there is a mix- 
ture of n — and n — 2 Landau levels (e.g. the state 
is of the form c + |2) + + c_|0)_ if the field is in the +z 
direction, where the coefficients will generally depend on 
temperature and field). Following Ref. [TSJ, we call these, 
respectively, the A and U phases. Writing out the Yang- 
Mills Lagrangian from Eq. ([I]) in terms of the fields w±, 
it can be shown that the coefficients in our theory are 
K\ = 1, K2 = and K 3 = —1, so it appears that we 
are on the boundary between the stable and unstable re- 
gions, and also-if we assume that the phase boundary 
remains unchanged for the holographic superconductor- 
on the boundary between the A and U phases described 
above. One can imagine changing the coefficients of these 
gradient terms by hand, thereby moving away from this 
critical point, but in that case we would no longer be 
dealing with pure Yang-Mills theory, which has the at- 
tractive feature of having no adjustable parameters. 

Of course, the extent to which these results are rele- 
vant for holographic superconductors ought to be ques- 
tioned. Since our ansatz leads to a Lagrangian in which 
there are no terms coupling gradients in the z direction 
to those in the xy plane, the criterion for stability should 
remain the same as in the non-holographic case. Also, 
by assuming that one order parameter component van- 
ishes and employing separation of variables, it is shown 
below that the A phase is a solution in the stable region 
of the phase diagram for the holographic case. Separa- 
tion of variables, however, cannot be used to obtain the 
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more complicated U phase described above. It is possi- 
ble that there is an analogue of this phase in the holo- 
graphic superconductor, but because the U phase in the 
non-holographic case is a complicated function involving 
multiple Landau levels, we expect that the same would be 
true in the holographic case. Characterizing such a state 
would most likely involve minimizing numerically the free 
energy of variational wavefunctions, and the problem of 
determining the exact form of such a state and compar- 
ing its energy with that of the A phase is an interesting 
problem left for future study. It is important to note, 
however, that the criterion given above for distinguish- 
ing the A and U phases will not necessarily hold in the 
holographic case, so it is possible that the A phase is in 
fact the unambiguous ground state for the particular La- 
grangian that defines our theory. We therefore proceed 
pragmatically, assuming stability and the existence of a 
state with iw_ = and following the original approach 
of Abrikosov [2U], which was also used in Refs. [T4l 121) . 
(If B < 0, the following discussion holds if we replace 
w+ — > W-.) 

Since there will not be rotational symmetry for the vor- 
tex lattice as there was for the isolated vortex, we once 
again take all fields to be functions of all three spatial 
variables, with A" = A^(x,y,z). In order to make nota- 
tion more transparent and to allow for easier comparison 
with the existing literature on superconductivity, we now 
switch to more conventional notation in which the su- 
perconducting order parameters are represented by two 
complex scalar fields. Letting 



171,2 = Al iV + iA 
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the Yang-Mills part of the action in Eq. ([TJ can be ex- 
pressed as 
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where we have assumed that we are in the A phase de- 
scribed above and sufficiently close to the upper critical 
field that we can set the second order parameter com- 
ponent 7/_ = 0. Here and in what follows, the gradient 
operator is defined as V = (d x ,d y ), and to simplify no- 
tation we have let Af — > A$. 

As was done in Refs. [IHHI], near the upper critical 
field, we can define e = (H c2 — H)/H C 2 and expand the 



fields: 

y, z) = $ (0) + e$ (1) + C(e 2 ) 
A x>y (x, y, z) = 4°) + eA^ + 0(e 2 ) (13) 
V+ (x, y, z) = e 1 ' 2 ^ + e 3 ^) + {e 5 ' 2 ). 

To leading order near H c2 , the electromagnetic fields are 
$ = fi(l - z/zo), A y = xH c2 and A% = 0. The higher 
order terms take into account the backreaction of the 
bosonic field on the electromagnetic fields. The magnetic 
field in the boundary theory is H = (d x A y — d y A x )\ z —o. 
Taking the equation of motion for rj + from Eq.( |12| and 
letting r/ + (x,y,z) = m + (x 1 z;p)e wv gives 



=d z (gd z m A 
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where we have neglected the term ~ |m + | 2 m + , since m + 
is small near H c2 . We note here that the vortex lattice 
solution for the p-w&ve superconductor was not obtained 
in Ref. [14] because the gauge field ansatz in that paper 
contained a complex phase factor, so that A^ ~ e lpy . 
In our work, it is the complex bosonic field r/ + which 
is given the phase factor, which, in the language of the 
original gauge fields, corresponds to A^ ~ g 2 pv t . This is 



necessary to obtain the term ~ {H c2 x + p) in Eq. (14 1 



The distinction is important because, as we described 
above, the gauge field in this theory is real, with the role 
of the real and imaginary parts of the usual Ginzburg- 
Landau order parameter being played here by the t 1 and 
t 2 directions in SU(2) space. 



Taking advantage of the linearity of Eq. (14), we can 



again use separation of variables, letting m + (x, z;p) = 
p(z)j(x;p). We then obtain the following eigenvalue 
equations: 
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where X = \/H c2 {x + p/H c2 ). The first of these is just 
a harmonic oscillator equation, which is solved by the 
Hermite polynomials: 



ln {x;p) = e- x ' 2 / 2 H n (X). 



(17) 



Here n = 0, 1, 2, . . . denotes the Landau energy level, and 
the corresponding eigenvalues are given by A„ = 2n + 1 . 
The Abrikosov vortex lattice is given by a superposition 
of the lowest energy [n = 0) solutions: 

m + (x,y,z) = p (z)^2c :j e lpjy 'jo(x;p j ), (18) 
i 

where the Cj are coefficients that determine the structure 
of the vortex lattice. As shown in Refs. [TH [T51 [2Tj . 
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FIG. 2. Upper critical magnetic field for the w + component 
of the superconducting order parameter. The lower left part 
of the phase diagram is the superconducting region, and the 
upper right is the normal state. Note the upward curvature, 
which is a direct consequence of the ^-dependence in Eq. ( 16 1. 



where integration by parts has been used, and a bar de- 
notes fields evaluated on shell. We assume that the scalar 
field n + has compact support in the (x, y) coordinates, 
so that the contributions from the boundaries x = const. 
and y = const, vanish when we integrate by parts. Fur- 
thermore, there is no boundary contribution from the 
horizon due to the regularity condition, and none from 
the AdS boundary at z = due to the boundary condi- 
tion n + (z — 0) = 0. 

Similarly, the on-shell conditions for the fields $ and Ai 
can be calculated from Eq. ( 20 ) and the result substituted 



back into the action, yielding the action with all fields 
evaluated on shell 



qYM 
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d 3 xA l d z A, 
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the upper critical field H c i can be calculated at a given 
temperature by finding the highest field at which Eq. 
( 16 1 has a non- vanishing solution, indicating the presence 
of a superconducting condensate. The resulting phase 
diagram is shown in Fig. [2] We alert the reader to the 
characteristic upward curvature of H c i- This curvature, 
stemming from Eq. (16 1, is intrinsic to our theory and is 



thus reflective of the effects of field-induced correlations 
captured within the holographic approach. 

We next investigate the nature of the vortex lattice 
solution near H C 2. It is well known that the free en- 
ergy in an ordinary (s-wave) Ginzburg-Landau theory is 
minimized when the vortex cores form a triangular lat- 
tice, and this has also been shown to be the case for a 
holographic s-wave superconductor [21]. Here we follow 
a similar approach to find the configuration that mini- 
mizes the free energy of the holographic p-wave super- 
conductor. While our analysis is complicated somewhat 
compared to the s-wave case by the quartic term and the 
many gradient terms that appear in Eq. (12), we shall 



find that-just as in the s-wave case-the triangular vor- 
tex lattice minimizes the free energy. 

Since all quantities are time-independent, the free en- 
ergy is given by SI = —Sq^ 1 / J dt, where Sqo 1 is the 
action evaluated with the fields n+, A XtV and <f> on shell. 



The equation of motion for n+ following from Eq. ( 12 1 is 



=d z (gd z n+) + ~(D X + iD y ){D x - iD y )n+ 

$ 2 1 , 

+ — V+ + (d x A y - d y A x )n+ - -\n+\ n+. 
9 * 
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Multiplying this equation by n* + and integrating over 
space, combining this result with Eq. ( 12 ) gives the ac- 
tion with ?7 + evaluated on shell: 



S YM (V + ) = 
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Here there is no boundary contribution from the $ terms 
upon integration by parts due to the compact support 
in the (x, y) coordinates and the boundary condition 
$(2 = 0) = 0. As discussed in Ref. [2T], the bound- 
ary contributions from the gauge fields at the horizon can 
also be ignored since they are independent of the field n + , 
and our interest here is in finding the configuration of 77+ 
that minimizes the free energy. The boundary term at 
z = for the field Ai does not in general vanish, however; 
and, according the AdS/CFT dictionary, the expectation 
value of the current operator in the boundary theory is 
proportional to d z Ai. 

The leading nonzero correction to the action is at 
C(e 2 ), so the free energy can be expressed as SI ~ 

51(0) +e 2^(2) with 
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For simplicity we no longer use the bar to denote fields 
that are on shell. In order to evaluate Eq. (22 1, we need 



to obtain an expression for A ^\ At 0(e), the equation 
of motion following from Eq. (12 I is 



[d z (gd z ) + 2V 2 ]A 
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We have also chosen the gauge condition diAi = and 
applied it in deriving Eq. ( 23 ) . It is important to note 



that Eq. ( 24 ) describes currents in the bulk theory, and 



is distinct from the boundary current operator alluded 
to above. Furthermore, since ?y+(z,x) ~ 7(x), which de- 
scribes the ground state of a simple harmonic oscillator, it 



is a simple matter to show [55] that j 
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The antisymmetric symbol satisfies 6j, = — e,-j and ei2 
1. 



The solution to Eq. ( 23 ) can be expressed as 
v4, (1) (z,x) = aj (x) 



/ dz 'dx! 'Gb(z, z ';x — x)j^ (z',x'), 



where a, (x) is the homogeneous part of the solution sat- 
isfying eijdidj — —H C 2, and the Greens function Gb sat- 
isfies 

[^ z ( 5 d z ) + V 2 ]G B (z,z';x-x') 

= -5(z - z')5(x - x') (26) 
Gb(z — 0, x) = lim g(z)GB(z, z'; x) = 0. 



We can now use Eq. ( 25 ) to evaluate the second term 
in Eq. p2|: 



dxAf ] d z A t 



J dz'dx'dxa l (x)9 z G B (z,z';x-x')jf ) (z',x') (27) 
J dz'dx!dxd z GB(z,z';x-x!)\ri ( +\z\x!)\ 2 , 



H C 2 
~2~ 



where in the second equality integration by parts and 
VGb(z,z';x — x') = —VGb(z,z';x — x') were used. 
From the boundary condition in Eq. ( 26 ) , the integral 
of Gs can be performed: 

rmin(z,z') j n 

dxG s (z,z';x) = / — — . (28) 
Jo 9\ z ) 



Using this result along with Eq. (|27|), we obtain the fol 
g expri 



lowing expression for the free energy: 
-1 
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As pointed out in Ref. [21], however, Eq. (29) is not 



our final result, since it depends on the normalization 
of T}+ . This ambiguity in normalization is resolved by 



considering nonlinearity. From Eq. (19 1, the equation of 



motion for r}+ is Crfl^ = 0, where we have defined the 
differential operator 

t =d z (gd z ) + \{D X + iD y )W(D x - iD y )W 



H 1- ti c2 . 



(30) 



Also following from Eq. ( 19 1 is the equation of motion 



for T)+\ which can be expressed as Crj^ — J, where 



J^ l -{A x +iA y )^{D x -iD y )^^ 
+ l (D x + W y )W \(A x -iA y )V v ( V 



( 9x A y -d v A x )U - *$W<Z>U + 
9 * 



(31) 



With this equation of motion, we obtain the following 
identity: 



= I d 3 x 



(l)* T (l)* A ( 2 ) 



= I d 3 



i]^*J - ( Ci] 

d 3 XT)W*J, 



(32) 



where integration by parts along with the boundary con- 
dition r] + (z — 0) = was used to obtain the sec- 
ond line. Again using integration by parts and j 



(i) 



\iijdj\ry + | 2 in Eq. (32), we obtain the condition 



d A x 



1 (33) 



In order to evaluate Eq. ( 33 1 we must determine the 
form of At 0(e), the equation of motion for $ is 



(34) 



( 5 ^ + v 2 )$« = | ?? i 1) | 2 $ (0) . 

This equation has solution 

<f (1) (z,x) = J dx! J dz'^^G t (z,z'; X - X ') 

where Gt(z, z'; x — x') is the Green function satisfying 

( ff <9 2 + V 2 )G 4 (z,z';x-x') = 
- 9 (z)(5(z-z')<5 (2) (x-x'). 



(36) 



In order to make further progress, we expand the Greens 
functions from Eqs. (26) and (36) in a basis of eigenfunc- 
tions: 



G t (z,z';x) = J2U^W)G2(x, A) 

A 

-g(z)d 2 z Uz) = A&(*) 
6(0) = = £ x (zo), 



(37) 



and 



G B (z,z';x) =^2xx(z)x{(z')G 2 (x,X) 
a 

~\dz [g{z)d z x\{z)\ = Axa(z) 
X a(^ = 0) = lim g(z) X '(z)=0. 

z->z 

where the Greens function G 2 satisfies 

(V 2 -A)G 2 (x,A) = -£(x). 



(38) 



(39) 



Due to Eqs. (|25[) and (35), the free energy in this 



theory takes a nonlocal form, as opposed to the usual, 



7 



non-holographic Ginzburg-Landau theory, which is com- 
pletely local. This is due to the fact that the Ginzburg- 
Landau theory is a low energy effective expansion, 
whereas the AdS theory presented here retains the 
physics from all energy scales [U]. To get a local ef- 
fective theory, we recognize that, in the long wavelength 
limit, G2(x, A) decays much more quickly than |?7^(x)| 2 , 
so we can approximate 



J dx'G 2 (x-x',A)| 7 (x')|' 



l7(x)| 2 
A ' 



(40) 



Using Eqs. (25 



form of the cone 
limit: 



and (35), we can now give an explicit 



ition Eq. (33) in the long wavelength 



= d 



1 2 ( 

4 ~~g 



r^V 2 / dVG B (2,z';x~x')hl 1) (2',x')| 2 (41) 
p 4 &(z)p 2 " 



|7| 4 -iW|7| 2 



To obtain the first equality we have again used integra- 
tion by parts and V'Gb(z, z 1 ; x — x') = -VGb(z, z'; x — 
x'), and the second equality gives the approximate form 
in the long wavelength limit, using Eqs. ( 37 )-( 39 ) . We 
have also defined 



rameter: 

Q _ I r 
V ~ V 

v 



2(H c2 - Hf (p 2 ) 2 (2p 4 - ap 2 ) 
q 2 K 2 /3 (p 4 — ap 2 ) 2 



(43) 



where (. . .) denotes spatial average, and the Abrikosov 
parameter is given by 



/3 



(I7I 4 ) 
(111 2 ) 2 



(44) 



Since the free energy density in Eq. ( 43 1 is negative, min- 
imizing the free energy corresponds to minimizing (3. It 
is well known that the vortex lattice distribution that 
minimizes f3 is the triangular vortex lattice, for which 
f3 = 1.159. This was also the result found for the holo- 
graphic s-wave superconductor |21j . 

In conclusion, we have shown the existence of a vortex 
solution in a holographic p-wave superconductor at low 
magnetic fields, as well as a vortex lattice solution near 
the upper critical magnetic field, H c2 . H c2 exhibits a 
characteristic upward curvature, intrinsic to our theory, 
which reflects the effects of field-induced correlations cap- 
tured by the holographic approach. The free energy was 
found to be minimized by the triangular vortex lattice. 
In the future it would be interesting to extend this theory 
to a BCS-like theory of fermions in AdS [33] , which would 
give insight into the possible types of p-wave pairing in 
holographic superconductors, as well as the tantalizing 
possibility of Majorana fermions, which are known to ex- 
ist as bound states in the vortex cores of chiral p-wave 
superconductors. 



a(z) 



4$(0) 



By combining Eqs. (41) and (29), we can now give an 



approximate, local expression for the free energy density 
that is independent of the normalization of the order pa- 
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